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ABSTRACT 

Using high-resolution Hubble Space Telescope data, we reexamine the fundamental 
properties (ages, masses and extinction values) of the rich star cluster population in 
the dwarf starburst galaxy NGC 5253. The gain in resolution compared to previous 
studies is of order a factor of two in both spatial dimensions, while our accessible 
wavelength range transcends previous studies by incorporation of both near-ultraviolet 
and near-infrared (IR) passbands. We apply spectral synthesis treatments based on 
two different simple stellar population model suites to our set of medium-, broad-band 
and Ha images to gain an improved physical understanding of the IR-excess flux found 
for a subset of young clusters (30 of 149). With the caveat that our models are based 
on fully sampled stellar mass functions, the NGC 5253 cluster population is dominated 
by a significant number of relatively low-mass (M c \ < a fewlO 4 M ) objects with ages 
ranging from a few x 10 6 to a few x 10 7 yr, which is in excellent agreement with 
the starburst age of the host galaxy. The IR-excess clusters are almost all found in 
this young age range and have masses of up to a few xlO 4 M Q . The IR excess in 
the relatively low-mass NGC 5253 clusters is most likely caused by a combination of 
stochastic sampling effects and colour variations due to the presence of either luminous 
red or pre-main-sequence stars. We also find a small number of intermediate-age (~ 1 
Gyr-old), ~ 10 5 M clusters, as well as up to a dozen massive, ~ 10 Gyr-old globular 
clusters. Their presence supports the notion that NGC 5253 is a very active galaxy 
that has undergone multiple episodes of star cluster formation. 

Key words: methods: data analysis - stars: luminosity function, mass function - 
galaxies: starburst - galaxies: star clusters: general - infrared: galaxies 



1 INTRODUCTION 
1.1 Choice of target 

NGC 5253, a metal-poor blue compact dwarf galaxy in the 
Centaurus A group (Karachentsev et al. 2007), hosts an ex- 
tremely young starburst. At the fortuitously close distance 
of 3.13 Mpc (Davidge 2007; 1 arcsec corresponds to 16 pc), 
the galaxy has been the target of a large number of obser- 
vational programmes spanning a range of wavelengths, from 
X-rays to the radio domain (e.g., Rieke, Lebofsky & Walker 
1988; Caldwell & Phillips 1989; Martin & Kennicutt 1995; 
Freedman et al. 2001; Cresci, Vanzi & Sauvage 2005; Martm- 
Hernandez, Schaerer & Sauvage 2005; Lopez-Sanchez et al. 
2007, 2011). 

* E-mail: grijs@pku.edu.cn 



The metallicity of NGC 5253, approximately 0.2-0.4 
Z Q (Kobulnicky & Skillman 1995; Kobulnicky et al. 1999; 
Lopez-Sanchez et al. 2011), is relatively constant across the 
galaxy, except for a few areas exhibiting higher-than-average 
nitrogen abundances (Walsh & Roy 1989; Kobulnicky et 
al. 1997; Lopez-Sanchez et al. 2007; Lopez-Sanchez & Es- 
teban 2010a,b; Monreal-Ibero et al. 2010) and a possible 
slight helium enrichment (Campbell, Terlevich & Melnick 
1986; Lopez-Sanchez et al. 2007). This low metallicity, com- 
bined with the galaxy's relatively small dimensions (~ 5 x 2 
arcmin 2 ~ 4.6 x 1.8 kpc 2 major/minor axis extents) and low 
mass (~ 2-5 x 10 9 M Q ; Bottinelli et al. 1972; Turner, Beck 
& Hurt 1997), render its starburst component an excellent 
analogue to galaxies in their early formation phases, while 
its proximity makes it ideal for studies of the many bright, 
young stellar clusters (YSCs; e.g., Gonzalez- Riestra, Rego 
& Zamorano 1987; Caldwell & Phillips 1989; Calzetti et al. 
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1997; Harris et al. 2004; Vanzi & Sauvage 2004; Cresci et al. 
2005; Martfn-Hernandez et al. 2005; Harbeck, Gallagher & 
Crnojevic 2012). 

In fact, based on ultraviolet (UV) Hubble Space Tele- 
scope (HST) /Space Telescope Imaging Spectrometer (STIS) 
spectroscopy of the YSCs and diffuse background light in the 
main body of the galaxy, Tremonti et al. (2001) suggested 
that star clusters may have been forming continuously, and 
subsequently dissolve on ~ 10-20 Myr timescales (based on 
arguments related to tidal effects worked out in detail in 
Kim, Morris & Lee 1999; a process since coined cluster 'in- 
fant mortality'), thus dispersing their stars into the field star 
population. 

1.2 Starburst characteristics 

The galaxy's Ha morphology (Martin 1998; Calzetti et al. 
2004; Meurer et al. 2006) shows that the inner core is cur- 
rently undergoing intense star formation and is, in fact, one 
of the youngest starbursts known (van den Bergh 1980; 
Moorwood & Glass 1982; Rieke et al. 1988; Caldwell & 
Phillips 1989; Calzetti et al. 1997; McQuinn et al. 2010a,b). 
The active starburst is confined to the inner 60 pc (Calzetti 
et al. 1997, 1999; Tremonti et al. 2001); age estimates for 
the starburst duration in this region are < 1 x 10 7 yr, and 
more likely a few xlO 6 yr, assuming a burst-like formation 
history for the diffuse stellar population in the disc (Walsh 
& Roy 1989; Rieke et al. 1988; Tremonti et al. 2001). These 
estimates are based on weak [Fe n] emission, strong Br7 
emission, a lack of significant (non-thermal) synchrotron 
emission from supernova remnants (Beck et al. 1996), the 
presence of Wolf-Rayet stars (Campbell et al. 1986; Kob- 
ulnicky et al. 1997; Schaerer et al. 1997; Lopez-Sanchez et 
al. 2007; Lopez-Sanchez & Esteban 2010a,b; Monreal-Ibero 
et al. 2010) and weak 2/xm emission, implying only a small 
contribution of giants and supergiants, which thus provides 
a rough upper age limit. 

In two very detailed studies of the galaxy's central star- 
burst region, Calzetti et al. (1997, 1999) showed that the 
UV emission is dominated by a 3-4 Myr old, relatively small 
star cluster, NGC 5253-4 (Meurer et al. 1995). On the other 
hand, a dust-enshrouded (Ay = 9-35 mag) very young cen- 
tral cluster (2-2.7 Myr old), NGC 5253-5, is responsible for 
most of the nuclear region's ionization (Tremonti et al. 2001; 
Cresci et al. 2005). At radio wavelengths, Turner et al. (1998, 
2000) derive ionization potentials driven by large numbers 
(200-1000) of O-type stars, which they interpret to imply 
that very large clusters are the preferred mode of star for- 
mation in the central regions of the galaxy. 

1.3 Current context 

It is thus clear that the galaxy's star and star cluster for- 
mation history is highly complex. Despite the plethora of 
existing, detailed studies, many of the salient details of the 
galaxy's evolution remain open-ended. This realization is 
perhaps mostly driven by the object's proximity, which al- 
lows us to probe the star-forming environment in signifi- 
cantly greater detail than in more distant starbursts. On 
the other hand, NGC 5253's proximity also offers us the key 
advantage of it having been observed repeatedly across the 
full observable wavelength range. 



Nevertheless, where the stellar component has been sub- 
ject to detailed scrutiny, most studies have either focussed 
on careful, mostly spectroscopic studies of a few objects at 
a time, or systematic exploration of the galaxy's full extent 
but based on relatively simple approaches. This is particu- 
larly so for its star cluster population. Previous statistical 
studies have relied on straightforward application of simple 
stellar population (SSP) analysis based on optical broad- 
band data (e.g., Harris et al. 2004; Cresci et al. 2005; Har- 
beck et al. 2012). However, in recent years, theory and com- 
putational approaches have made significant advances. 

In this paper, we go significantly beyond previous stud- 
ies of the NGC 5253 star cluster population in a number of 
complementary ways. First, although we still use primarily 
(medium- and) broad-band imaging observations, we base 
our results on the highest-achievable spatial resolution from 
the near-UV (NUV) to the near-infrared (NIR) regimes. The 
gain in resolution compared to previous studies is of order a 
factor of two in both spatial dimensions, while our accessible 
wavelength range transcends previous studies by incorpora- 
tion of the crucial NUV and NIR passbands (which will allow 
us to at least partially break the age— met allicity and age- 
extinction degeneracies). Where relevant, we also refer to 
Ha emission-line results to further reduce the uncertainties. 

Second, we focus specifically on one aspect of extra- 
galactic star cluster populations where we take full advan- 
tage of the observational gains, i.e., the importance and ori- 
gin of a red or IR excess in the spectral-energy distributions 
(SEDs) of a relatively low-mass cluster population such as 
that in NGC 5253. To address this aspect, we do not need to 
have access to a statistically complete sample of star clusters 
down to the lowest detection limits enabled by the observa- 
tions' worst signal-to- noise ratio (S/N). 

This paper is organized as follows. In Section 2 we dis- 
cuss the observational data which form the basis of our anal- 
ysis, as well as our data reduction approach. The analysis 
methods adopted are outlined in Section 3, while we discuss 
the robustness of the results and perform a comparison in 
Section 4. We offer a preliminary discussion of the effects 
of stochasticity in the stellar mass function separately, in 
Section 5, and provide an overview of the physical state of 
the NGC 5253 cluster population in Section 6. Finally, we 
summarize and conclude the paper in Section 7. 

This is the first article in a series of three. In Paper II 
(P. Anders et al., in prep.), we will present a detailed theo- 
retical investigation of the effects of stochastic sampling on 
cluster age and mass determinations based on adopting fully 
sampled stellar mass functions using our own SSP models 
and analysis approach (cf. Section 3). In Paper III (R. de 
Grijs et al., in prep.), we will apply these insights to the data 
set analysed here, to explore in detail to which extent our 
results are affected by biases owing to stochastic sampling 
of the clusters' mass functions. The present paper is meant 
to establish the baseline for this body of work. 

2 HUBBLE SPACE TELESCOPE 

OBSERVATIONS AND DATA REDUCTION 

2.1 Data and basic reduction 

We mined the HST Data Archive for available medium- 
/broad-band and Ha observations of NGC 5253 obtained 
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with its main optical and NIR cameras, i.e., the Advanced 
Camera for Surveys (ACS), the Wide-Field and Plane- 
tary Camera-2 (WFPC2) and the Near-Infrared Camera 
and Multi-Object Spectrometer (NICMOS). Since we had 
selected our target galaxy based on the availability of a 
large number of high-resolution HST data sets spanning the 
longest possible wavelength range, we could make a careful 
selection of the best possible data sets to reach our science 
goals. Table 1 includes details of the final set of observations 
used, which were selected to achieve the best possible com- 
bination of spatial resolution (aimed at reducing contami- 
nation by neighbouring sources), S/N (for which we initially 
used exposure time as a proxy) and extensive, continuous 
wavelength coverage. 

2.2 Object selection 

We used the ACS/HRC data in the F330W, F435W, F550M 
and F814W filters as basis for the final orientation of all ob- 
servations. This implies that we rotated, scaled and aligned 
all other images to the position, scale and orientation of 
the ACS/HRC data using the iraf/stsdas tasks 1 rotate, 
magnify, and imalign, respectively, using a selection of 
conspicuous sources common to each frame as guides. We 
checked that our image manipulation did not affect the re- 
liability of the integrated flux measurements of our sample 
of candidate star clusters. The final image size of the set 
of aligned images is 1169 x 1138 ACS/HRC pixels, 2 corre- 
sponding to 32.73 x 28.45 arcsec 2 (~ 500 x 430 pc 2 ). 

The standard deviations (c 8 k y ) of the number of counts 
in seemingly empty sections in all images were established 
to ascertain a 'sky' background count for each filter. Multi- 
ples of this background count were used as thresholds above 
which the numbers of 'real' sources in both the F550M and 
F814W ACS/HRC filters were calculated, using the idl 3 
find task. We did not force our detection routine to con- 
strain the sources' roundness or sharpness, in order to be as 
inclusive as possible. The most suitable thresholds for source 
inclusion, at 4<j s k y (see below for justification), were 0.7 
and 1.5 counts s _1 (270 and 248 source detections) for the 
F550M and F814W filters, respectively. In all passbands, the 
number of detections initially decreases rapidly with increas- 
ing threshold value. This is an indication that our 'source' 
detections are noise-dominated in the low-threshold regime. 
Where the rapid decline slows down to a more moderate 
rate, our detections become dominated by 'real' objects (ei- 
ther stars, clusters or cluster candidates, or real intensity 
variations in the background field: for an example, see fig. 3 
of Barker et al. 2008). In the remainder of this paper we will 

1 The Image Reduction and Analysis Facility (iraf) is dis- 
tributed by the National Optical Astronomy Observatories, which 
is operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the U.S. Na- 
tional Science Foundation. STSDAS, the Space Telescope Science 
Data Analysis System, contains tasks complementary to the ex- 
isting IRAF tasks. We used Version 3.6 (November 2006) for the 
data reduction performed in this paper. 

2 The ACS/HRC pixels are rectangular and subtend 0.028x0.025 
arcsec 2 each. 

3 The Interactive Data Language (idl) is licensed by Research 
Systems Inc., of Boulder, CO, USA. 



only consider the objects in the 'source-dominated' domain, 
with fluxes in excess of the relevant 4<r s k y threshold. 

We subsequently employed a cross-identification pro- 
cedure to determine how many sources coincided with 
intensity peaks within 1.4 pixels of each other in the 
F550M®F814W HRC filter combination (i.e., allowing for 
1-pixel mismatches in both spatial directions). In our next 
step, we applied a Gauss-fitting routine in idl to each can- 
didate cluster. We also determined the best-fitting Gaus- 
sian profiles for some 10-12 isolated, star-like sources in all 
filters (where the number of objects available for this pur- 
pose depended on the S/N of the observations and the typ- 
ical stellar spectra) and thus obtained observational stellar 
Gaussian <tg (size) values for all filters. The Gaussian stellar 
widths in the F550M HRC band (the filter we used for our 
size selection) were found to be best represented by a Gaus- 
sian profile of era — 1.1 pixels. Any candidate cluster with 
fG,F550M < 1.1 pixels was considered most likely to be a star 
(or - if much smaller than 1 pixel - an artefact either of the 
detector or caused by cosmic rays) and was consequently 
discarded. We applied a conservative clipping criterion so 
as not to reject some marginally extended sources. These 
steps resulted in rejection of 88 objects from the original 
F550M source list and led to a working sample containing 
182 extended cluster candidates. 

Their size distribution is shown in Fig. 1, where we also 
show the linear cluster sizes after deconvolution of the ob- 
served ct g and the intrinsic size of the point-spread function 
(PSF). Our (candidate) cluster sample is characterized by 
a peak at sizes close to that of the PSF, while the most 
extended cluster candidate has <tg = 3.4 pixels (equivalent 
to 1.36 pc after PSF deconvolution). Our sample is insen- 
sitive to stellar associations because of the selection crite- 
ria adopted. We will discuss the implications of the cluster 
size distribution (as well as the insets in Fig. 1) in Section 
6.4. In a final step, we confirmed that our procedure is ro- 
bust with respect to misidentifications of genuine clusters as 
point sources (stars): only a few (<5) per cent of the ob- 
jects that were rejected based on the Gaussian size estimates 
were brighter than the tip of the red giant branch (TRGB) 
magnitude at the distance of NGC 5253, where we used the 
most recent calibration based on observations of lu Centauri, 
m trgb = _4.04±0.12 mag (Bellazzini et al. 2001, 2004). 

We note that although realistic cluster luminosity pro- 
files may well deviate from the simple Gaussian profile 
adopted here, its consistent and systematic application to 
our observations allows us to differentiate accurately be- 
tween objects of different sizes, irrespective of their true pro- 
files, provided that the intrinsic profiles of the individual ob- 
jects do not differ too much from object to object (although 
one should realize that stochastic sampling effects may in- 
validate this assumption; see below). In Fig. 2 we show the 
distribution of our candidate cluster sample across the face 
of the NGC 5253 common field in the ACS/HRC F550M 
filter. 4 The main differences between our current sample se- 
lection and that of Cresci et al. (2005) - the most recent 
detailed statistical study of the NGC 5253 cluster system 
- are owing to the higher-resolution images we used as the 



4 For a movie including all high-resolution images and object 
identifications, see http://kiaa.pku.edu.cn/~grijs/ngc5253.mov. 
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Table 1. HST observations used for the analysis in this paper. 



Filter 


Camera 


Chip Exposures (s) 


Total cxp. 


Proposal 


PI 


Photometric corr. 








time (s) 


ID 




factor (unit) 



F330W 


ACS 


HRC 


3 x 449 


F435W 


ACS 


HRC 


4 x 150 


F547M 


WFPC2 


WF 


2 x 200, 2 x 600 


F550M 


ACS 


HRC 


4 x 200 


F555W 


ACS 


WFC 


4 x 600 


F606W 


WFPC2 


PC 


500 


F658N 


ACS 


HRC 


240 


F814W 


ACS 


HRC 


4 x 92 


F110W 


NICMOS 


NIC2 


4 x 23.97 


F160W 


NICMOS 


NIC2 


4 x 23.97 


F222M 


NICMOS 


NIC2 


4 x 63.95 



1347 


10609 


Vacca 


1.0 


(cthrc) 


600 


10609 


Vacca 


1.1 


(CHRC) 


1600 


6524 


Calzctti 


1.6 


(o"wpc) 


800 


10609 


Vacca 


1.0 


(<7hrc) 


2400 


10765 


Zczas 


1.0 


(ctwfc) 


500 


5479 


Malkan 


2.0 


(f HRC) 


240 


10609 


Vacca 


1.0 


(°"HRC) 


368 


10609 


Vacca 


1.25 (fTHRC ) 


95.86 


7219 


Scovillc 


1.2 


(ctwfc) 


95.86 


7219 


Scoville 


1.1 


(<T W fc) 


255.80 


7219 


Scovillc 


1.15 (o- W fc) 



Note: The three-digit numbers in the filter names represent their effective wavelengths in nanometres; the final digits indicate the filter 
type: N = narrow, M = medium, W = wide. 
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Figure 1. Gaussian size distribution of the extended star cluster 
candidates in our sample (as fraction of the total sample). The 
arrow indicates our minimum size cut-off, corresponding to the oq 
of point sources. The top axis shows the linear cluster sizes after 
deconvolution of the observed ctq an d the intrinsic size of the PSF. 
The insets show the size distribution (top) in logarithmic units, 
combined with the best-fitting power-law function (dashed line, 
defined by the proportionality given) and (bottom) as a function 
of F547M magnitude (see Section 6.4 for discussion). The vertical 
dotted line indicates the ACS/HRC PSF size. 



basis for our object identification, thus reducing crowding 
effects. As a result, we find slightly fewer cluster candidates 
in the galaxy's main body; some of the sources identified by 
Cresci et al. (2005) as 'clusters' are unresolved (and likely 
stars) in our images (with ctg < 2.65 pc, corresponding to a 
PSF-deconvolved size of ffG.dec < 2.4 pc). 



2.3 Photometry 

Our custom-written idl aperture photometry task uses 
source radii and sky annuli individualized for each cluster 
candidate. The adopted radii and annuli were based on 
Gaussian stellar size measurements (ctg) obtained in the 
F550M filter. We applied a correction factor as a function of 
filter/instrument to correct for different intrinsic PSF sizes 
(see the final column in Table 1), based on determination 
of the best-fitting Gaussian profiles for the 10-12 isolated, 
star-like sources in all filters we used for our size cuts in Sec- 
tion 2.2. We used the F550M HRC and F555W ACS/Wide 
Field Camera (WFC) images as our basis for the determi- 
nation of source and sky annuli because of their high reso- 
lution and good S/N: in images with HRC (or similar, e.g., 
WFPC2/PC) resolution, we used a source aperture radius of 

4 x [(Tf55om x HRC correction factor] , and 4 and 6 x [ctfssom x 
HRC correction factor] for the inner and outer sky annuli 
(in pixels), respectively. For the other images, characterized 
by poorer resolution, the relevant numbers were 3, 3.5 and 

5 x [ctf555w x WFC correction factor] , respectively, again ex- 
pressed in pixel units (the WFC chips are characterized by 
a pixel size of 0.049 arcsec). The units in the final column 
of Table 1 corresponding to cthrc and <twfc refer to ofssom 
and OF555W, respectively. 

This choice of apertures and their scaling with the ob- 
jects' measured sizes was based on extensive tests in which 
we inspected the stellar radial profiles, to identify where 
the object profiles generally disappear into the background 
noise. Specifically, we carefully checked the profiles of all 
candidate clusters, for all filters and detectors, to verify that 
our source radii were chosen conservatively and sufficiently 
far out so as not to exclude any genuine source flux with re- 
spect to the general background level (i.e., well in excess of 
the radii where the radial profiles disappear into the back- 
ground noise) . We also checked that the background annuli 
were chosen appropriately and not dominated by neighbour- 
ing sources. Finally, we checked that the background level 
was generally flat (within the photometric uncertainties and 
the statistical noise) so that small differences among the 
radii used to define the background annuli as a function of 
filter would not cause artificial offsets in flux, and hence lead 
to colour differences. 
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Figure 2. NGC 5253 field of view in common of all observations and source selection discussed in this paper, overlaid on the ACS/HRC 
F550M image. Red squares: Age t ~ 10 10 yr; blue triangles: t ~ 10 9 yr (for a discussion of the clusters' spatial distribution as a function 
of age, see Section 6.1). 



Fig. 3 shows a typical example of the level of photomet- 
ric accuracy attainable with individualized apertures com- 
pared to the use of fixed apertures. For this example, we 
employed source apertures of 8 pixels (in the detector's na- 
tive system) and background annuli between 8 and 12 pix- 
els for our ACS-based fixed-aperture photometry. For the 
F547M WFPC2/WF and the NICMOS images, we used 
equivalent radii of 12, 12 and 15 pixels, again in the respec- 
tive detectors' native systems. The photometric data points 
in the different panels only include those objects for which 
the individualized apertures were smaller than or equal to 
the sizes of our fixed apertures, so that a direct one-to-one 
comparison is possible. The mean differences in the sense 
(wimdiv. ap . — infixed ap .), as well as the standard deviations 
of these difference distributions are included in each of the 
panels for the relevant instrument/filter combination. It is 
clear that the correspondence between our individualized 
aperture photometry and that obtained from using fixed 
apertures (chosen so that most sources would have most 
of their flux recorded) is tight. The scatter increases for the 
lower-resolution and lower-S/N images, although the mean 
difference remains small, particularly within the associated 
uncertainties. 

For the size cut-offs (l.kr G ) in the HRC F550M and the 
WFC F555W filters, the annuli adopted in this study corre- 
spond to deconvolved linear sizes of 2.35, 2.35 and 3.50 pc, 
and 2.45, 3.85 and 4.10 pc, respectively. The resulting instru- 
mental magnitudes of the 182 'clusters' for which our tasks 
could obtain reliable photometry (with photometric uncer- 
tainties of < 0.3 mag in at least 7 filters) were adjusted for 
their zero-point offsets (determined from the PHOTFLAM 
and PHOTZPT header keywords of the HST observations, 
thus resulting in photometry in the ST mag system). This 
sample forms the basis of our analysis in Sections 4 and 5. 
We will make the full photometric data tables of all 182 



cluster candidates available in Paper III, upon completion 
of our detailed investigation of the effects of stochasticity in 
the clusters' stellar mass functions. 

Note that because of the extended nature of our objects 
(which was, after all, a key selection criterion), we could not 
apply PSF photometry or simple aperture corrections (but 
see, e.g., Harris et al. 2004; Annibali et al. 2011); the latter 
would require uncertain assumptions about, e.g., the un- 
derlying source profile (cf. Anders, Gieles & de Grijs 2006) 
as well as the degree of mass segregation within the clus- 
ters. The rationale for such an approach is to attempt to 
obtain the full flux of a given object. Particularly for the 
lower-luminosity, lower-mass objects, the effects of stochas- 
tic sampling would likely introduce additional uncertainties 
if adopting a generic underlying radial profile. On the other 
hand, adopting an aperture correction to large radii based 
on the profile and colours in a restricted region relatively 
close to a cluster's core may give rise to additional colour- 
and age-dependent uncertainties in the derived total fluxes 
owing to the effects of stellar mass segregation. Instead, we 
opted to use a hybrid method: our aperture sizes were vari- 
able, i.e., based on the measured cluster sizes and designed 
such that our output photometry included the full flux con- 
tribution from all objects. 5 

In addition, our integrated SEDs are dominated by 
the highest-luminosity (inner) regions of the clusters (we 
checked that any contribution from potentially missed flux 
in the outer regions is negligible), so that the resulting pa- 
rameters do not strongly depend on our aperture photom- 

5 In a sense, whether or not one should (or even could) use aper- 
ture corrections to recover the fluxes of extended objects is a 
philosophical issue. In this paper, we decided to adopt the lessons 
learnt from the results of Anders et al. (2006) and choose source 
apertures based on the objects' size measurements. 
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28 26 24 22 20 IB 28 26 24 22 20 18 28 26 24 22 20 18 

Magnitude (individualized apertures) (ST mag) 

Figure 3. Comparison of aperture photometry based on our in- 
dividualized apertures and that from the use of fixed apertures. 
For this example, we employed source apertures of 8 pixels (in the 
detector's native system) and background annuli between 8 and 
12 pixels for our ACS-bascd fixed-aperture photometry. For the 
F547M WFPC2/WF and the NICMOS images, we used equiva- 
lent radii of 12, 12 and 15 pixels, again in the respective detec- 
tors' native systems. The photometric data points in the different 
panels only include those objects for which the individualized 
apertures were smaller than or equal to the sizes of our fixed 
apertures, so that a direct one-to-one comparison is possible. The 
mean differences in the sense (mj nd j v ap — m,fl xct j ap ), as well 
as the standard deviations of these difference distributions are 
included in each of the panels for the relevant instrument /filter 
combination (also indicated). The dashed lines indicate the loci 
of equality. 



etry approach (see also Section 4.4). Finally, the majority 
of our sample clusters are sufficiently isolated that our ap- 
proach is not severely affected by the effects of crowding; this 
is supported by inspection of the growth curves of the dozen 
sample clusters in the most crowded regions. Ultimately, the 
main advantage of our hybrid approach is that it is simple 
and requires few, if any, assumptions that may introduce 
additional uncertainties. 

Fig. 4 compares the resulting photometry with that 
published by Harris et al. (2004). The latter authors pub- 
lished ST magnitudes for 33 objects they identified as clus- 
ters, 25 of which are also included in our cluster sample. We 
note that the overall agreement between both data sets is 
reasonable, although in half a dozen cases there are signif- 
icant offsets. The latter occur particularly for the brighter 
objects in the F547M filter. 

Harris et al. (2004) discuss the best approach to ob- 
taining reliable photometry in the crowded environment of 
NGC 5253, which is additionally affected by significant back- 
ground variations. They conclude that PSF-like photometry 
is inappropriate under these conditions, and that aperture 
corrections depend sensitively on the underlying luminosity 
profile adopted. They settled on achieving a balance between 




20 19 18 17 

P814..H04 (mag) [WFPC2/WF] 



Figure 4. Comparison of our photometry using a hybrid ap- 
proach (indicated by the subscript 'new') with that of Harris et 
al. (2004), based on aperture corrections appropriate for point 
sources (referred to as H04), both in the ST mag system. 



the systematic errors imposed by small apertures and the 
random errors inherent to adopting large apertures by se- 
lecting an intermediate aperture size of 5 pixels (radius). 
They then attempt to correct for the clusters' extended 
profiles by adopting a 0.1 mag, 5-pixel aperture correction 
appropriate for point sources, but they caution that it is 
likely that there remains a systematic error in the photom- 
etry because of the clusters' non-spherical, extended mor- 
phologies. In this context, the offsets in the F547M pho- 
tometry between our determinations and those of Harris et 
al. (2004) can be understood fairly easily. For the brigher, 
extended sources in particular, the 5-pixel apertures they 
adopted are too small to encompass the full source fluxes, 
thus leading to underestimated luminosities (fainter magni- 
tudes) on average. This is evident from inspection of Fig. 
4, particularly for objects with mF5547M < 20 mag (in both 
sets of photometric measurements). In addition, two bright 
objects (mF5547M,H04 < 18 mag) have brighter magnitudes 
in the data set of Harris et al. (2004) than determined by 
us. These objects are found in the most crowded regions 
of the galaxy, where their extended profiles overlap with 
neighbouring objects. On balance, we prefer to use our own 
newly determined photometry, because this attempts to in- 
clude and correct for the effects of non-sphericity and spatial 
extent. We also point out that a comparison of our photom- 
etry based on different apertures (shown in Fig. 3) indicates 
that our measurements are internally consistent. 

The scatter in the F814W photometry is somewhat 
larger than that for the F547M data, which is most likely 
caused by the difference in resolution between the cam- 
eras used in the presence of fairly crowded fields. Lower- 
resolution observations, such as the WFPC2 data of Harris 
et al. (2004), are more significantly affected by close neigh- 
bours ('blending') than the higher-resolution ACS data we 
use in this paper. In addition, there is a small systematic 
offset towards brighter magnitudes for our photometry com- 
pared to the Harris et al. (2004) data set, which is most 
likely caused by their underestimated aperture corrections 
for point sources. 

In Section 3 we analyse our cluster photometry based 
on the objects' colours. In essence, we use a poor man's ap- 
proach to cluster spectroscopy based on the clusters' broad- 
band SEDs. 
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3 SED ANALYSIS 

We base our SED analysis on two of the current-best, most 
up-to-date SSP model suites available, specifically the Yg- 
gdrasil models (Section 3.2) and the GALEV models (Section 
3.3). Both sets of models are based on different, modern 
isochrones and independently developed spectral synthesis 
codes, which have all been extensively validated. Although 
there are currently at least a handful of suitable SSP model 
suites to choose from, the two model sets we adopted are 
most appropriate in the context of the present paper, since 
a number of our co-authors are actively involved in their 
development. This offers the key advantage that we fully 
understand the basic underlying assumptions and essential 
physical inputs, which we summarize in detail in this section. 

3.1 General model assumptions 

For a one-to-one comparison between the observed and 
model SEDs, we assumed that star clusters are coeval 
(single-burst) stellar populations characterized by a Kroupa 
(2001) IMF. This IMF is characterized by a lower cut-off 
mass of 0.15 M Q , whereas the upper-mass cut-off ranges 
between ~ 66 and ~ 76 M , and is determined by the 
mass coverage of the Padova isochrones (i.e., Marigo et al. 
2008) for a given metallicity, e.g., the upper mass cut-off for 
Z = 0.004 (0.008) is 68.07 (70.35) M . The theoretical stel- 
lar libraries are from Lejeune et al. (1997, 1998) for a broad 
range of metallicities. 

The current versions of both sets of models include the 
important thermally pulsing asymptotic-giant-branch (TP- 
AGB) evolutionary phase. At ages from ~ 100 Myr to ~ 1 
Gyr, TP-AGB stars account for 25 to 40 per cent of the bolo- 
metric light, and for 50 to 60 per cent of the A"-band emis- 
sion of SSPs (see Chariot 1996; Schulz et al. 2002). The effect 
of including the TP-AGB phase results in redder colours for 
SSPs with ages between ~ 10 8 and 10 9 yr, with the strongest 
effect (up to > 1 mag) in (V — K) for solar metallicity, and 
in (V — I) for Z > 0.5 Z Q . Shorter- wavelength colours and 
lower- metallicity SSPs are less affected. 

The models predict the contributions to the combined 
SED of both the stars and the ambient, photo-ionized gas 
owing to the presence of young massive stars. The photo- 
ionized gas contributes both to the continuum and to specific 
emission lines, thus directly affecting the integrated broad- 
band fluxes of YSCs. Anders & Fritze-v. Alvensleben (2003) 
and Adamo et al. (2010b) showed that nebular emission non- 
negligibly affects cluster SEDs during the first 10-15 Myr 
of their evolution. Although this contribution reduces after 
6 Myr, at least for shorter (bluer) wavelengths, it remains 
significant in the NIR regime up to ~ 20 Myr (depending 
on metallicity). 

3.2 Yggdrasil models and IR-excess analysis 

Model-specific details. The primary approach we adopted for 
our analysis in this paper uses the most up-to-date Yggdrasil 
spectral synthesis models (Zackrisson et al. 2011). The lat- 
ter are based on Padova-AGB SSPs from the Starburst99 
spectral synthesis code (Leitherer et al. 1999; Vazquez & 
Leitherer 2005). For our application to the NGC 5253 clus- 
ter system, we adopted metallicities of Z = 0.004 and 0.008 



(where Z = 0.019). Adopting the Calzetti et al. (2000) 
attenuation law, which is characterized by total-to-selective 
absorption Rv = 4.05 and most appropriate for use in star- 
burst environments, we created a grid of models including 

0. < E(B - V) < 3.0 mag for each age step. 

For the gaseous contribution, our models employed 
standard Hn region values as input parameters to CLOUDY 
(version 90.05; Ferland et al. 1998), i.e., a hydrogen gas den- 
sity of 10 2 cm~ 3 , a filling factor of 0.01 and a covering fac- 
tor of unity (no leakage of ionizing photons, so that - for 
instance - all Lyman continuum photons contribute to ion- 
izing the gas). In addition, we assumed that the metallicities 
of the stars and gas in a given cluster are identical. 

Parameter determination. We used a least-squares fit 
to estimate, for each cluster, the best-fitting model (ex- 
pressed in magnitudes; see, e.g., Adamo et al. 2010a) and 
calculated the associated Q value (or x probability func- 
tion: models with Q>0.1 are good, while Q > 0.001 is still 
marginally acceptable; Press et al. 1992), as well as the re- 
duced xl = xV^> where v = N — m — lis the number of 
degrees of freedom given N data points and m fit param- 
eters. Here, N refers to the number of filters and m = 2, 

1. e., we obtain best fits of the age and extinction (m = 2) 
and a scaled mass of each cluster for fixed metallicity. The 
resulting 68 per cent confidence limits on the determination 
of the cluster parameters were based on the prescription of 
Lampton, Margon & Bowyer (1976). Adamo et al. (2010a,b, 
201 la, b) used the spectral synthesis models of Marigo et al. 
(2008) to validate their results. 

Adamo et al. (2010a) discovered that several young clus- 
ters in the luminous blue compact galaxy (BCG) Haro 11 
have a flux excess in the /, H and K s bands. Inclusion of 
those filters in the fit resulted in significant residuals and a 
clear overestimation of the cluster age. Given the starburst 
nature of NGC 5253, we expected that at least some of the 
galaxy's YSCs might be affected by a similar red excess. 

Therefore, we produced three series of y 2 minimiza- 
tions, i.e., one including all available data points from the 
NUV regime (F330W) to the F222M filter, a second exclud- 
ing F110W, F160W and F222M (i.e., from the NUV to the / 
band), and a third which also excluded the / band (leaving 
the range up to the F606W filter's effective wavelength). Our 
initial analysis showed that the latter filter selection did not 
produce any improvement to the fit (in the \ 2 sense) and 
this filter selection leaves too few constraints on the result- 
ing cluster parameters because of a number of degeneracies 
affecting the data, so that we discarded those results. 

We checked whether we could identify J-band excesses, 
but exclusion of the /-band filter did not produce any im- 
provement of the fit in the x 2 sense (as expected if the fit 
including the /-band photometry is correct). In this paper, 
we will thus only discuss the NGC 5253 cluster fit results 
based on application of our approach to (i) the full wave- 
length range and (ii) the range excluding the NIR filters. 

3.3 GALEV models and AnalySED approach 

Model- specific details. Our second, complementary approach 
uses the GALEV SSP models (Kotulla et al. 2009; and refer- 
ences therein, as well as subsequent unpublished updates). 
They cover ages between 4 x 10 6 and 12.6 x 10 9 yr, with 
an age resolution of Alog(i yr -1 ) = 0.05, and are based on 
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the latest set of stellar isochrones of the Padova group for 
metallicities of 0.0001 < Z < 0.03, tabulated as 15 discrete 
metallicities. Our model grid was completed by inclusion 
of extinction effects (Calzetti et al. 2000), with E(B - V) 
spanning the range from 0.0 to 2.0 mag, and adopting a 
resolution of 0.05 mag. 

Parameter determination. In recent years, we developed 
a sophisticated tool for star cluster analysis based on broad- 
band SEDs, AnalySED, which we tested extensively both 
internally (de Grijs et al. 2003a,b; Anders et al. 2004) and 
externally (de Grijs et al. 2005), using both theoretical and 
observed star cluster SEDs. 

Note that in this paper we base our main analysis on the 
Yggdrasil models and analysis approach discussed in Section 
3.2. We apply the ANALYSED approach to the exact same 
filter set as used in the IR-excess approach. This ensures that 
the results from both approaches are directly comparable. 

For a given cluster, the shape of the observed SED 
(again, expressed in magnitudes), including the relevant ob- 
servational uncertainties, was compared with the shapes of 
model SEDs (as a function of age, metallicity and foreground 
extinction) for the passband combination considered. Each 
model SED (and its associated physical parameters) was 
assigned a probability based on the \ 2 value of this compar- 
ison. The model data set with the highest probability (i.e., 
the lowest \ 2 value) was adopted as the most representative 
set of cluster parameters. Models with decreasing probabili- 
ties were summed up until a 68.26 per cent total probability 
(la confidence interval) was reached, to estimate the uncer- 
tainties in the best-fitting model (cf. Anders et al. 2004) . 



4 COMPARISON OF PARAMETER 
DETERMINATIONS 

In the remainder of this paper, we will use our Yggdrasil 
models and associated analysis, and the resulting fundamen- 
tal cluster parameters, as the basis for the ensuing discus- 
sion. In this section, we explain in detail how we determined 
which clusters are affected by an IR excess and pursue an 
approach in which we determine the basic cluster parame- 
ters. In Sections 4.1 and 4.2, we use the Yggdrasil models 
and analysis as our basis; in Section 4.3, we compare these 
results with those from application of the GALEV models and 
the AnalySED approach. 

4.1 Determination of IR excesses 

Based on application of our least-squares fitting routine to 
the full wavelength coverage available (i.e., all 10 filters with 
central wavelengths from 330 nm to 2.22 pm), we obtained 
best fits with sufficiently small xt an d Q > 0.001 for 149 of 
the 182 objects in our sample with photometric uncertain- 
ties < 0.3 mag in at least 7 filters. In the remainder of this 
paper, we will base our analysis on this reduced data set of 
149 clusters and discard the 33 objects for which we could 
not achieve a satisfactory fit (we will briefly return to these 
objects below). A number of clusters returned best fits with 
xl 3> 1 and Q < 10 -5 . We inspected the residuals, i.e., the 
differences between the observed SED and the model fit, in 
each of the N filters used for the fit (for all clusters, indepen- 
dent of the goodness-of-fit parameters) and confirmed that 



30 objects were characterized by an IR excess with respect 
to the best-fitting model that included all NUV and optical 
filters up to and including F814W. The remaining 119 ob- 
jects exhibited SEDs which could be fitted with 'standard' 
SSP models over the full wavelength range. 

Fig. 5 shows comparisons of the resulting cluster param- 
eters. On the horizontal axes, we plot the ages, extinction 
values and masses based on fits to the full F330W— S-F222M 
data set (10 filters), adopting a metallicity of Z = 0.004 
henceforth. On the vertical axes, we plot the results for the 
F330W-s>F814W 7-filter set. The 30 clusters with apprecia- 
ble IR excess are shown as the objects contained in the (red) 
circles. Overall, there is a reasonable one-to-one correlation 
between the results from the two different data sets. As we 
already pointed out in de Grijs et al. (2005), the mass de- 
terminations are clearly better constrained than the age or 
extinction determinations. The latter are very sensitive to 
small changes (including those due to photometric problems) 
in the SEDs. 

4.2 Systematic effects for all versus IR-excess 
sources 

Both the age and the mass comparisons for the sources ex- 
hibiting an IR excess show that exclusion of the NIR points 
results in smaller values. For the mass determinations, this is 
expected at least to some extent, because masses are deter- 
mined based on a scaling of the entire SED; excess intensity 
in a number of filters will hence artificially increase the scal- 
ing upwards. In addition, redder SEDs will result in older 
ages, and hence the resulting mass estimates will be based 
on incorrect mass-to-light ratios, because the latter are age- 
dependent. For the age determinations, the sources affected 
by an IR excess are artificially given older ages (and older 
clusters are assigned higher masses for the same luminosi- 
ties), assuming fully sampled cluster mass functions. 

For the full sample of 149 objects, we found that the 
exclusion of the NIR filters in our Yggdrasil-based analysis 
resulted in 112 objects that exhibited an age difference of 
Alogi[yr] < 0.05 dex. Eight objects were characterized by 
age differences of 0.05 < A log t [yr] < 0.30 dex, 18 objects re- 
sulted in 0.30 < Alogi[yr] < 0.50 dex and for seven objects 
Alogi[yr] > 1.0 dex. The equivalent numbers for the differ- 
ences in mass are 92 objects with A log M c i[Mq] < 0.05 dex, 
27 exhibiting 0.05 < A log M c i[M ] < 0.30 dex, 16 for which 
0.30 < AlogM c i[M ] < 0.50 dex and only three showing a 
logarithmic age difference in excess of 1.0 dex. 

We also checked whether our choice of metallicity 
may have affected the results through the well-known age- 
metallicity degeneracy. We refitted our entire data set, again 
using the two different filter combinations, but now adopt- 
ing Z = 0.008 for all clusters. The ages resulting from the 
Z = 0.008 assumption were very similar to but marginally 
younger than those based on Z = 0.004 (not shown). This 
effect translated in a small to negligible offset that did not 
appear to be a function of age. Thus, we conclude that our 
results are negligibly affected by the age-metallicity degen- 
eracy (but see Section 4.3), and that adoption of a single 
metallicity for all clusters does not appreciably affect our 
analysis. 

In Fig. 6 we compare the residuals of the fits obtained 
from single-mindedly fitting the entire available wavelength 
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All filters, [F330W, F222M] 



Figure 5. Comparison, using the models discussed in Section 3.2, 
of the fundamental cluster parameters as a function of wavelength 
coverage. X axis: all 10 niters; y axis: all 7 NUV and optical filters, 
up to and including F814W. The objects identified with (red) 
circles are characterized by an IR excess (see text for discussion) . 
Note that our standard SED analysis returns similarly young ages 
for many of the cluster candidates (cf. Fig. 7, top left-hand panel). 
Therefore, the density of data points in the top left-hand panel 
appears lower than it is in reality. 



range (shaded histograms) and those resulting after exclu- 
sion of the NIR passbands (open histograms drawn using 
solid lines), again for the clusters affected by an IR excess 
only. It is clear that the residuals in the latter case are more 
symmetrically distributed around zero. Note that the most 
significant effect is seen in the residuals of the F814W filter. 
These residuals improve upon exclusion of the NIR data be- 
cause this allows us to obtain a 'proper' match (in the \ 2 
sense) of all data points in the optical regime. Because of 
the flux excess in the NIR filters, in the fits including all 10 
filters the observed F814W flux was consistently below the 
best model prediction. 

Fig. 7 is a collection of the distributions of the basic 
cluster parameters for the full data set of 149 clusters. We 
show both the distributions initially obtained by fitting the 
entire available wavelength range (dotted histograms) and 
those resulting from NUV to optical filters only (solid lines) . 
As in Fig. 5, we see that this filter selection leads to younger 
ages and lower masses, while the overall distribution of ex- 
tinction values remains statistically unchanged. Note that 
the large majority of our sample clusters, particularly those 
exhibiting an IR excess (see also Fig. 5), have masses be- 
tween a few x 10 2 and 10 5 Mq (for confirmation of this mass 
range, see e.g., Harris et al. 2004; Cresci et al. 2005). These 
are most likely significantly affected by stochastic sampling 
of their stellar mass functions. Table 2 includes our final 
set of derived cluster parameters based on the Yggdrasil 
SSP models and associated analysis (assuming fully sam- 
pled mass functions). 
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Figure 6. Fit residuals as a function of filter for the clusters 
affected by a significant IR excess. The shaded histograms are 
the residuals based on fits to the full available wavelength range; 
the solid histograms show the results after exclusion of the filters 
affected by the IR excess. 



We also explored whether our results are affected sig- 
nificantly by the age-extinction degeneracy, which often 
causes artificial offsets if one has access to only broad-band 
SEDs covering smaller wavelength ranges. The optical-NIR 
colour-colour diagram of Fig. 8 clearly shows that the clus- 
ters in which we detected an IR excess are all located red- 
ward of the model SSP (blue solid line). This, therefore, 
serves as an independent check of the reality of our IR-excess 
determinations (see also fig. 2 of Cresci et al. 2005). 

The large majority of IR-excess sources are located to- 
wards blue optical colours. If they were affected by signif- 
icant extinction, we would have expected them to be dis- 
placed well along the direction of the reddening vector. How- 
ever, since we have access to up to 10 data points spanning 
the entire NUV/optical/NIR range, this helps our efforts 
to break the age-extinction degeneracy quite significantly: 
for instance, the evolution in age and the extinction vector 
are not oriented parallel to each other in the colour-colour 
diagram of Fig. 8, at least for young ages (they would be 
in optical-only colour-colour diagrams). We thus conclude 
that the young IR-excess sources are, by and large, negli- 
gibly affected by the age-extinction degeneracy and, more 
importantly, that the IR excess in these sources is most likely 
not owing to the effects of extinction. For reference, we also 
included the 33 objects for which we could not obtain satis- 
factory fits using fully sampled SSP models in Fig. 8 (small 
blue crosses). Their distribution in colour-colour space fol- 
lows that of the remaining 149 objects which are the main 
focus of this paper; a fraction also appears to exhibit an IR 
excess with respect to the standard SSP models. 

4.3 One-to-one comparisons: dependence on the 
SSP model suite adopted 

Finally, we compare the resulting masses and ages from our 
Yggdrasil models and analysis with those obtained based on 
the galev+ AnalySED approach: see Fig. 9 and Table 3. 
The top row of Fig. 9 shows the one-to-one comparisons of 
the clusters' ages and masses (the resulting extinction values 
cannot be determined as reliably). While the overall trends 
show reasonable agreement, particularly for the masses (cf. 
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Table 2. NGC 5253 cluster coordinates and derived parameters assuming fully sampled Yggdrasil SSP models." 



ID 


R.A. (J2000) 
(dog) 


Doc. (J2000) 
(dog) 


°"G,F550M 
(pixels) 


^filters' 


log(t) 

[yr] 


Uncertainties 
(positive) (negative) 


log(Af c i) 
[M ] 


Uncertainties 
(positive) (negative) 





204.9766060 


-31.6426963 


1.154 
















1 


204.9775826 


-31.6427056 


1.165 


10 


7.18 


0.67 


0.40 


3.13 


0.26 


0.61 


2 


204.9810049 


-31.6366573 


1.194 


10 


7.15 


0.03 


0.37 


3.04 


0.13 


0.46 


3 


204.9810436 


-31.6368153 


1.158 


10 


9.79 


0.33 


0.48 


4.67 


0.21 


0.39 


4 


204.9792423 


-31.6411761 


1.154 


10 


6.78 


0.03 


0.03 


2.71 


0.10 


0.03 



a The full data table is available associated with the online version of this article. This stub is meant to provide guidance regarding 
form and content of the full table. 6 Number of filters used to obtain the derived parameters. 




Figure 7. Distributions of the fundamental NGC 5253 cluster 
parameters. The dotted histograms with dashed shading repre- 
sent the parameters based on fits to the full wavelength range 
from the F330W to the F222M filter; the histograms drawn us- 
ing solid lines are the resulting distributions after removal of the 
NIR filters from the SED analysis. The bottom right-hand panel 
shows the distribution of our sample clusters in the age-extinction 
plane (see Section 6.4). The objects identified with red circles are 
characterized by an IR excess. 



de Grijs et al. 2005) as well as for most of the IR-excess 
sources, there are a few discrepancies (indicated by the dot- 
ted circles) that require further investigation. 

Fig. 10 displays the SEDs of the four clusters for which 
our Yggdrasil-based analysis and the ANALYSED approach 
obtained the most different age estimates [At (yr) > 2a]. 
The latter returned an age for cluster 107 (indicated by the 
small dotted circle in the top left-hand panel of Fig. 9) of 
2.0±0.8 Gyr and an extinction value of E(B — V) = 0.0 mag, 
while our Yggdrasil-based analysis resulted in estimated age 
and extinction values of 6.0 ±0.7 Myr and 0.60 mag, respec- 
tively, for similar-quality fits. Clusters 100, 130 and 167 (con- 
tained in the large dotted circle in the same panel) are also 
clearly discrepant because of the age-extinction degeneracy: 
the (age, extinction) pairs for each of these clusters are, re- 
spectively, (4.5 Myr, 1.05 mag) versus (11.2 Gyr, 0.25 mag), 



Figure 8. Colour— colour diagram in the ST mag system. Blue 
heavy solid versus magenta heavy dashed line: GALEV versus Yg- 
gdrasil SSP models for Z = 0.004, spanning an age range from 3 
Myr (blue extremeties) to ~ 12 Gyr (red extremeties of the mod- 
els). Green points: Theoretical stellar loci (Marigo et al. 2008); 
small blue crosses: Objects for which we could not obtain a sat- 
isfactory fit using fully sampled SSP models. The objects identi- 
fied with red circles are characterized by an IR excess. The black 
dashed line delineates the region in colour-colour space where star 
clusters and individual stars may be found (taking into account 
realistic photometric uncertainties). 



(5.0 Myr, 1.10 mag) versus (11.2 Gyr, 0.25 mag), and (4.5 
Myr, 0.90 mag) versus (12.6 Gyr, 0.05 mag), where the first 
and second (age, extinction) pairs refer to the results from 
the ANALYSED and Yggdrasil-based analyses, respectively. 
Additionally, differences between both sets of models and 
the details of the fitting routines may contribute to a lesser 
extent. The significantly different age estimates for these 
four clusters are associated with similarly varying mass-to- 
light ratios, which in turn also lead to significant differences 
in the resulting mass estimates. This effect is seen in the top 
right-hand panel of Fig. 9, where the objects of relevance are 
indicated by the large dotted circle. 

Despite these differences for a small number of clusters, 
the main features in the overall age and mass distributions 
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Table 3. NGC 5253 cluster parameters derived assuming fully sampled GALEV SSP models. 2 



11 



ID log(t) 



Uncertainties 



log(M cl ) 



Uncertainties 



[yr] 


(positive) 


(negative) 


[M ] 


(positive) 


(negative) 


6.65 


0.05 


0.05 


2.62 


0.11 


0.03 


6.65 


0.55 


0.05 


3.24 


0.18 


0.32 


6.60 


0.15 


0.05 


2.70 


0.32 


0.03 


7.00 


0.35 


0.05 


3.73 


0.48 


0.10 



a The full data table is available associated with the online version of this article. This stub is meant to provide guidance regarding 
form and content of the full table. 
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Figure 10. SEDs (absolute magnitudes in the ST mag system) of clusters 100, 107, 130 and 167 (contained in the dotted circles in 
Fig. 9), showing the effects of the age-extinction degeneracy (see text). The model SEDs are all based on GALEV models, adopting the 
best-fitting ages and extinction values derived from the ANALYSED and Yggdrasil-bascd ('IR-excess') approaches. 



(bottom row) are fairly well reproduced by both independent 
approaches. Assuming a simple underlying skewed Gaussian 
distribution, the peak in the age distribution resulting from 
our AnalySED analysis (including all objects) occurs at 
(log(t yr -1 )) = 7.11, with ff(iogt) = 0.80. For comparison, 
the Yggdrasil-based analysis returned (log(t yr -1 )) = 7.10, 
with cr/iogt) = 0.94. For the mass distributions, we find 
(log(M c i/MQ)) = 3.38 versus 3.37, and o-( logMd ) = 0.76 ver- 
sus 0.77, respectively. This excellent correspondence is much 
better than the feature-reproduction robustness reported in 
de Grijs et al. (2005), where we quoted reproducibilities of 
A(log(i yr -1 )) < 0.35 and A(log(M c i/M )) < 0.14. 

If we compare the age and mass determinations for 



the individual clusters, a similar result emerges. The 
mean differences between the Yggdrasil-based analysis and 
the AnalySED results is (Alog(t yr -1 )) = 0.041 and 
(Alog(M c i/M0)) = 0.100, respectively, with spreads (Gaus- 
sian cr's) of 0.900 and 0.517 dex, respectively. 



4.4 To which extent does Ha photometry make a 
difference? 

So far, we have focussed on the power of (medium- and) 
broad-band SED analysis based on extensive wavelength 
coverage. This is the type of observational data one might 
expect to obtain most easily for any extragalactic star clus- 
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Figure 9. (top) Comparison of the individual cluster ages and 
masses based on application of our two different SSP-matching 
approaches. The objects identified with (red) circles are charac- 
terized by an IR excess. The dotted circles indicate the handful 
of sources which are most likely significantly affected by the age- 
extinction degeneracy (see Section 4.3). (bottom) Age and mass 
distributions based on application of the ANALYSED approach. 

ter system, and hence exploration of the limitations of such 
a data set is relevant. For our analysis of the specific star 
cluster sample of NGC 5253, we also obtained narrow-band 
F658N observations from the HST Data Archive, covering 
the Ha emission line. 6 For most of the NGC 5253 clusters, 
numerous authors, including Calzetti et al. (1997) and Cresci 
et al. (2005), provide support for the young ages resulting 
from their fits by referring to the Ha fluxes found associated 
with these objects. Ha fluxes have also been found useful for 
analyses of numerous other YSC systems (for recent refer- 
ences, see e.g., Chandar et al. 2011; Whitmore et al. 2011; 
Fouesneau et al. 2012; and references therein), although the 
present analysis represents the first full integration of Ha 
photometry in a simultaneous multi-passband (N > 4 fil- 
ters) approach. 

In this section, we add the F658N/Ha photometry to 
our SEDs for confirmation and reanalysis of the properties 
of the NGC 5253 clusters using the Yggdrasil models. We 
obtained two sets of F658N photometry for our star clus- 
ter sample. First, we adopted the same apertures as for the 
nearest continuum filter, F550M. However, Ha emission may 
be more extended than that of the associated continuum 
sources. As such, we determined the Gaussian ctg's ° n the 
F658N image itself and used these size estimates to deter- 
mine a second set of F658N magnitudes (using the source, 

6 Note that NGC 5253 was also observed through a narrow- 
band filter covering the H/3 line. Calzetti et al. (1997) already 
extensively explored the galaxy's morphology and dust proper- 
ties based on these observations, so that we will not repeat that 
analysis here. 



Figure 11. (a) F658N photometry as a function of aperture size 
based on the intrinsic Gaussian ctq's in the F550M and F658N 
images using the individualized apertures defined in Section 2.3. 
Red data points: objects younger than 10 Myr; dotted line: lo- 
cus of equality, (b) Cluster photometry in the F550M and F658N 
filters, using aperture sizes based on oq fssom • Dotted line: lo- 
cus of equality corrected for the difference in filter transmission 
windows, (c) and (d) Mass and age comparisons, respectively, ob- 
tained using the photometry based on the differently sized aper- 
tures from panel (a). We have left out the relevant error bars 
for presentational clarity. All panels include the data or derived 
properties for all 149 clusters in our final sample. 

inner and outer sky radii expressed in units of ctg as defined 
in Section 2.3). A comparison of Gaussian sizes (not shown) 
indicates that the Ha emission coincident with our sample 
clusters exhibits a much broader range than the associated 
continuum emission in the F550M filter. 

Fig. 11a shows the resulting F658N photometry as a 
function of aperture size used. Although a reasonable one- 
to-one correlation between the F658N magnitudes resulting 
from using aperture sizes based on the F550M and F658N 
<tg's is apparent, there is a clear offset to brighter magni- 
tudes if the F658N-based size criterion is used, particularly 
for objects returned as younger than 10 Myr, based on the 
Yggdrasil-based broad-band analysis (indicated in red) . This 
underscores the need to account for differences in source 
sizes as a function of passband. 

Panel (b) shows the distribution of our sample clusters 
in the F550M versus F658N magnitude plane (the latter 
includes both continuum and line flux, if present). Again, 
the red data points are the objects younger than 10 Myr 
and the dotted line is the locus of equality, corrected for the 
difference in filter transmission windows. This panel shows 
that all sources that are more than twice their photometric 
la uncertainties from the dotted line are young (< 10 Myr). 
Note that we used the same apertures here for both axes, 
so that we can properly compare the presence of any Ha 
excess. 

Figs 11c and d compare the masses and ages, respec- 
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Figure 12. (top row) One-to-one comparisons for the full sample 
of 149 clusters of the ages, masses and extinction values based on 
using our medium- and broad-band photometry only ('BB'), and 
that based on also including the F658N data, (bottom row) Dot- 
ted histograms: distributions based on the medium- and broad- 
band data only; solid histograms: based on additional inclusion 
of the F658N fluxes. 

tively, obtained using the photometry based on the differ- 
ently sized apertures shown in panel (a). We have left out 
the error bars for presentational clarity (see, e.g., Fig. 5 for 
guidance) . These panels show that the impact of the differ- 
ent aperture sizes is very small, even at ages where signifi- 
cant Ha emission may be expected. 

In Fig. 12 (top row) we compare the ages, masses and 
extinction values based on using our (medium- and) broad- 
band photometry only ('BB'), and that based on also includ- 
ing the F658N data. The bottom row shows the equivalent 
histograms of these parameters. The dotted histograms are 
the distributions based on the (medium- and) broad-band 
data only; the solid histograms are based on additional inclu- 
sion of the F658N fluxes. The differences, in any of the pan- 
els, are very small and statistically negligible. This is exem- 
plified by a quantitative comparison of the histograms' char- 
acteristics. The peak and Gaussian width of the age distri- 
bution resulting from our analysis including the F658N pho- 
tometry are (log(i yr -1 )) = 7.13 and <?(\ og t) = 0.94 dex, re- 
spectively, compared with 7.10 and 0.94 dex for the analysis 
based on the (medium- and) broad-band images only. Simi- 
larly, the mean and Gaussian width of the equivalent mass 
distributions are (log(M c i/Mo)) = 3.36 and o"(i og M cl ) = 0.77 
dex, respectively, versus 3.37 and 0.76 dex. We thus con- 
clude that the Ha photometry fully supports our age and 
mass determinations of the NGC 5253 cluster sample based 
on medium- and broad-band photometry alone. 



5 INITIAL ASSESSMENT OF THE EFFECTS 
OF STOCHASTICITY 

Stochastic sampling of the stellar mass function, particu- 
larly of the massive stars in relatively low-mass ( < a few 
xlO 4 M Q ) - and therefore relatively poorly populated - 
clusters, may significantly affect the determination of fun- 
damental cluster parameters based on broad-band SEDs 
(e.g., Cervino, Luridiana & Castander 2000; Cervino et al. 
2002; Cervino & Luridiana 2004, 2006; Barker et al. 2008; 
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Mafz Apcllaniz 2009; Popescu & Hanson 2010; Fouesneau 
& Lancon 2010; Fouesneau et al. 2012; Popescu, Hanson & 
Elmegreen 2012). In general, stochastic sampling of cluster 
stellar mass functions leads to broad, asymmetric colour dis- 
tributions compared with the equivalent distributions asso- 
ciated with fully sampled SSP models and, as a consequence, 
important local biases are introduced in cluster age deter- 
minations for masses M c \ < a few xlO 4 Mq. The resulting 
cluster masses are usually less affected by stochastic sam- 
pling than the corresponding age estimates; they tend to 
scatter around the equivalent mass determinations derived 
using the traditional approach. 

Inspection of the derived cluster masses in Figs 7 and 
9 clearly indicates that the majority of our sample clusters 
may be affected by stochastic sampling of their mass func- 
tions (cf. the mass ranges quoted in Harris et al. 2004; Cresci 
et al. 2005) . This observation is further compounded by the 
results of our age determinations (based on the assumption 
of fully sampled SSPs!), which imply that a large fraction of 
the NGC 5253 clusters have ages near 10 7 yr (as confirmed 
and tightly constrained by the Ha analysis in Section 4.4). 

This corresponds to stellar population ages when lumi- 
nous red stars are present. Red-supergiant stars (RSGs; with 
stellar masses up to ~ 40 M Q ), in particular, start to appear 
around this time in realistic stellar populations, while at the 
same time the massive, hot, ionizing stars that dominate at 
young ages may cause significant variations in the nebular 
emission from stochastically sampled clusters (e.g., Cervino 
et al. 2000; Fouesneau et al. 2012). These stars are sig- 
nificantly brighter than their main-sequence counterparts, 
but not proportionally more massive. The absence or pres- 
ence of a small number of such stars can significantly affect 
the observed integrated colours of a cluster. For instance, 
stochastic effects are expected to affect the colour and mag- 
nitude of a 10 Myr-old, 10 5 Mq cluster by AV = 0.08 and 
A(V - /) = 0.07 mag (Dolphin & Kennicutt 2002). In ad- 
dition, these sources are relatively short-lived, and so the 
luminosity of the clusters might be expected to vary signif- 
icantly on short timescales (e.g., Kiss et al. 2006; Szczygiel 
et al. 2010; and references therein). This is of particular in- 
terest and relevance in the present context because our ob- 
servational data were taken over a period spanning a total 
of ~ 11 years between the first and final images included in 
this study (proposal IDs 5479 and 10609, respectively). 

We can obtain a better (initial) handle on the signif- 
icance of the stochastic sampling effects associated with 
our cluster sample by examining Fig. 8 in detail. First, 
we note that most of our clusters are young and relatively 
unaffected by extinction. For Ay = mag, it appears 
that most of the clusters can be explained by a young age 
( < 10 8 yr). However, in this case neither clusters redder than 
(m.F435w — wif550m) = 0.4 mag, nor those redder than the 
(mF8i4W — fiFi60w) values defined by the SSP models can 
be explained. This implies that at least some amount of ex- 
tinction is required to match the distribution of the clusters 
in colour-colour space. In addition, clusters that are red- 
der in (mF8i4W — f7iFi60w) than the envelope of the SSP 
models with variable extinction are most likely explained by 
the presence of individual stars. This includes the majority 
of objects we identified as having an IR excess. (Note that 
we cannot exclude the possibility that, in a small number 
of cases, an additional photon source that is unaccounted 
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for may also give rise to their violation of the observational 
boundaries for genuine stars and SSPs.) 

As such, we conclude that the IR excesses observed in 
30 of the NGC 5253 clusters are consistent with the effects of 
stochastic sampling. Similarly, the additional clusters found 
in the same part of the colour-colour diagram are likely 
also affected by stochastic sampling if we adopt the same 
reasoning. This includes 30-40 per cent of the total sample 
discussed in this paper. 

A first investigation focussing on the accuracy of derived 
cluster ages obtained along similar lines as done in this paper 
was published by Silva- Villa & Larsen (2011). We are in the 
process of implementing direct stochastic modelling using 
the GALEV SSPs as our basis. A preliminary comparison of 
our models with the results of Silva- Villa & Larsen (2011) in- 
dicates good agreement (E. Silva- Villa, priv. commun.). We 
note, however, that we are introducing important improve- 
ments with respect to their work. In particular, our obser- 
vations of the NGC 5253 cluster system cover a much larger 
wavelength range, while we are also investigating the effects 
of extinction as a free parameter in our fits. A much more in- 
depth investigation of the effects of stochastic sampling on 
the ages and masses determined based on broad-band pho- 
tometry will be published in Paper II. Once we have fully 
validated our stochastic modelling approach, we will under- 
take a detailed study of the effects of stochastic sampling on 
the derivation of the NGC 5253 fundamental cluster param- 
eters (Paper III). This will have important consequences for 
our understanding of the system, given that many clusters 
have masses in the regime that is dominated by stochastic 
sampling effects. The current cluster parameters will then 
provide a very useful comparison data set. 



6 PHYSICAL STATE OF THE NGC 5253 
CLUSTER POPULATION 

6.1 The cluster age distribution 

Fig. 13 shows the distribution of our final, well-fitted NGC 
5253 star cluster sample (of 149 objects) in the diagnostic 
age-mass plane, as well as the age and mass distributions 
resulting from our IR-excess analysis. The cluster population 
is dominated by a significant number of relatively low-mass 
(M c i <C 10 5 Mq) objects. Assuming fully sampled stellar 
mass functions, these objects have ages from a few xlO 6 
to a few x 10 7 yr. This age range is in excellent agreement 
with the starburst age of the host galaxy, < 10 Myr, and is 
confirmed by our age redeterminations based on the clusters' 
Ha fluxes. 

Our analysis has also revealed the presence of a small 
number of intermediate-age (~ 1 Gyr-old), ~ 10 5 Mq clus- 
ters, as well as up to a dozen old clusters resembling true 
globular clusters (GCs), in relation to both their ages (~ 10 
Gyr) and masses (> 10 5 Mq). We emphasize that our dis- 
covery of a small number of old GCs is secure. First, in our 
analysis we use a low metallicity of Z = 0.004 (0.2Zq). It is 
expected that the oldest stellar population components in 
a given galaxy are characterized by the lowest metallicity, 
since they were formed at a time when there was little recy- 
cled gas in the interstellar medium. (This argument is valid 
under the assumption that these clusters formed in NGC 
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Figure 13. Age— mass distribution of the NGC 5253 clusters. The 
objects identified with (red) circles are characterized by an IR ex- 
cess. The solid line in the age-mass plane represents our average, 
galaxy-wide detection limit at My — —5.4 mag (V ~ 22.1 mag), 
based on the SSP models adopted. For reasons of presentational 
clarity, we do not include error bars on the individual data points. 



5253 and were not subsequently accreted; see below.) Our 
age estimates at the current low metallicity are already al- 
most as old as the oldest ages in our SSP models, whereas 
at such old ages SSP models are quite insensitive to changes 
caused by metallicity differences. We confirmed that adopt- 
ing an even lower metallicity will not result in a different 
conclusion. 

In Fig. 13, we have also indicated the expected detec- 
tion (50 per cent completeness) limit based on standard SSP 
analysis (solid line); the limit is a steeply increasing func- 
tion of cluster mass with increasing age. Although we did 
not set out to select a statistically complete sample of star 
clusters, it is satisfying to note that our sample clusters obey 
a similar, galaxy-wide average detection limit across the en- 
tire age range. We focus specifically on resolved yet compact 
clusters (for which it is straightforward to obtain the inte- 
grated photometry and a generic completeness limit) rather 
than on more dispersed stellar associations (cf. Annibali et 
al. 2011). This also enables a better comparison with pre- 
vious studies of this galaxy, many of which were based on 
lower-resolution images than used in the present study. 

The selection of the final sample shown in Fig. 13 is gov- 
erned by our observations with the lowest S/N. Depending 
on the nature of the clusters (i.e., with or without an IR ex- 
cess), the images with the lowest S/N are the F330W and the 
F222M observations, respectively. The complex behaviour 
of the selection limit resulting from the data reduction and 
analysis steps, combined with the galaxy's highly variable 
background, renders determination of the 'true' detection 
limit less than trivial. In addition, the effects of stochastic- 
ity will also play a role in the (non-)detection of the clus- 
ters closest to the detection limit (e.g., Silva- Villa & Larsen 
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2010, 2011). Keeping in mind both of these issues, so that it 
is unrealistic to assign a single detection limit to the entire 
galaxy, the distribution of the data points in Fig. 13 implies 
that the average, galaxy-wide detection limit is My ~ —5.4 
mag (V ~ 22.1 mag), but we note that local variations, 
combined with the effects of stochasticity, could be up to 
0.3 mag. 

Given the number and concentration of old GCs above 
our detection limit in the NGC 5253 data, we speculate that 
these are merely the surviving high-mass, high-luminosity 
clusters of an initially much larger population of clusters 
that formed around the time of galaxy formation, many of 
which may have been disrupted or naturally faded to be- 
low the detection limit because of stellar evolution. The 
young clusters are distributed in two narrow 'chimneys' near 
log(i yr -1 ) = 6.7 and 7.2, which are caused by the features 
in SSP models around that age (cf. the 'wiggles' in the SSP 
model used to indicate the selection limit in Fig. 13), but are 
not necessarily associated with stochastic sampling of their 
mass functions. In reality, we would expect these clusters to 
be spread out more evenly across this young age range. 

The clusters in both chimneys are distributed ran- 
domly throughout the galaxy, similarly for objects with 
both log(t yr" 1 ) < 7.0 and 7.0 < log(t yr" 1 ) < 8.0. The 
YSCs characterized by a clear IR excess are confined to the 
galaxy's main body; the intermediate-age clusters in NGC 
5253, with ages \og(t yr -1 ) ~ 9, are found in the galaxy's 
outer region (cf. Harbeck et al. 2012), towards the south. 

Finally, the complement of > 10 Gyr-old clusters are 
predominantly located in the southeast of the main body of 
NGC 5253 (cf. Fig. 2). Since our observational data reach 
similar depths across the galaxy (except for the central re- 
gion on the East featuring a wide dust lane), we would 
have expected to detect such old massive clusters across 
the galaxy with a similar level of confidence. As such, we 
conclude that this spatial bias is likely real. 

If these clusters were moving on chaotic, freely tumbling 
orbits through the galaxy, one would have expected them to 
be well mixed on these long timescales. However, Kobul- 
nicky & Skillman (1995) showed, based on observed Hi gas 
dynamics, that the main body of the galaxy appears to be 
rotating as a solid body along the galaxy's major axis. Un- 
der these circumstances, clusters that formed at one end of 
the gaseous central bar or disc, may have remained confined 
to a small area since the time of their formation. Alterna- 
tively, the old clusters could have been accreted externally: 
to explain the galaxy's peculiar velocity field, Kobulnicky 
& Skillman (1995) also suggested that accretion of a low- 
metallicity, gas-rich companion galaxy on a higly inclined 
orbit is highly plausible (but see Davidge 2007). This in- 
terpretation is supported by the unusual weakness of the 
galaxy's CO emission, the morphology of the CO gas and 
its unusual kinematic properties (Turner et al. 1997). 

6.2 Red- and infrared-excess sources 

We also note that almost all clusters affected by an 
(infra-)red excess, are found among the YSC population 
in clusters with masses <10 4 M . Adamo et al. (2010a,b, 
2011a, b) identified two possible reasons for the occurrence 
of an (infra-)red excess in more massive clusters in their 
sample BCG galaxies. For clusters younger than approxi- 



mately 6 Myr, they suggested that the IR excess could be 
due to diffuse, hot dust in which these YSCs may still be 
embedded and/or a large fraction of young stellar objects 
and pre-main-sequence (PMS) stars. In the optical I band 
(F814W), the shallower red excess observed in such YSCs 
may be owing to a dust photolumincscence bump in their 
SEDs at 0.7-0.9 /im, often seen in the presence of very strong 
UV fields. Older clusters, with ages > 10 Myr, may exhibit 
an IR excess, which these authors attribute to the uncer- 
tain treatment of RSG stars in SSP models or the effect 
of stochasticity in the stellar mass function. In NGC 5253, 
we have probed to much lower masses than Adamo et al. 
(2010a, b, 2011a,b) could reach in their more distant BCG 
sample. 

Most of our IR-excess clusters are found in the younger 
'chimney', characterized by an age of ~ 5 x 10 6 yr. This 
is close to the lower age limit of our SSP models and also 
close to the age when YSCs are expected to become optically 
visible (i.e., after dispersion of at least some of the embed- 
ded dust associated with their formation process), so that a 
number of these clusters may have even younger ages. Based 
on the detailed studies referenced above, these clusters are 
most likely affected by a combination of stochastic sampling 
effects and the presence of PMS stars. The IR-excess clusters 
in the older chimney, at log(t yr -1 ) = 7.1-7.2 are also of rel- 
atively low mass and, hence, we expect them to be affected 
by stochastic sampling effects, while at this age the stochas- 
tic colour variations owing to RSG stars will also become 
apparent. 

Only two of our IR-excess clusters allow us to offer a less 
ambiguous explanation. The highest-mass cluster affected 
by an IR excess has a mass of almost 10 5 M . Although it 
has been assigned an age of ~ 4x 10 6 yr, stochastic sampling 
should only play a minor role, so that we are confident that 
the excess in this case is caused by the presence of a hot, 
diffuse dust cocoon or perhaps young stellar objects. The 
oldest cluster affected by an IR excess has an age of almost 
10 s yr. Although it has a low mass of ~ 6 x 10 3 M , so that 
stochastic sampling effects must duly be taken into account, 
this is also the time when (relatively short-lived) AGB stars 
appear in realistic stellar populations, playing a similar role 
as RSG stars at ages near 10 7 yr, both in terms of the num- 
bers present at a given time and their short lifetimes. 

6.3 Comparison with previous studies based on 
fully sampled stellar mass functions 

Cresci et al. (2005) published the most recent statistical 
study of the NGC 5253 cluster population (see also Har- 
ris et al. 2004). They based their results on HST F547M 
and F814W broad-band observations and ground-based K s - 
band observations with the Very Large Telescope, comple- 
mented with HST-based Ha data. Qualitatively, our results 
are similar to theirs, although there are a number of impor- 
tant differences. They divide their final cluster sample of 115 
sources into a 'young' subsample of 51 objects which were 
detected in the Ha filter, and an 'old' subsample. 

For the 'young' subsample, they derive ages spanning 
the range from 3 to 19 Myr. This corresponds to the bulk of 
the clusters in our two young chimneys in Fig. 9. To derive 
the cluster masses, they assumed an average extinction of 
Ay = 0.98 mag - which appears reasonable for most clusters 
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based on our extinction estimates in Fig. 7 and assuming a 
Calzetti et al. (2000) extinction law - and an average age 
of 8 Myr; for the 'old' clusters they assumed an average 
age of 20 Myr (cf. Harris et al. 2004). Above their 50 per 
cent completeness limit, they find a steeply declining number 
of clusters with increasing mass, with a maximum mass of 
~ 10 6 Mq. Our mass function (Fig. 9) samples significantly 
lower-mass clusters. 

Contrary to their analysis, however, we believe that we 
do not have sufficient numbers of clusters to comment on 
the reality (or otherwise) of a 'peak' in the mass function. 
Their result may be an artefact of the steeply increasing in- 
completeness limit with increasing age which could mask the 
true underlying mass distribution and produce an artificial 
turnover in the observed mass distribution (cf. Fig. 13). 

Prior to their work, Calzetti et al. (1997) used opti- 
cal HST/WFPC2 broad- and narrow-band observations in 
a careful attempt to characterize the star cluster proper- 
ties in the core of NGC 5253. Tremonti et al. (2001) used 
HST/STIS UV spectroscopy of four of the brightest clusters 
in the galaxy's central region to derive their fundamental 
properties. Of the four clusters examined in detail in both 
the Calzetti et al. (1997) and Tremonti et al. (2001) studies, 
we managed to produce satisfactory results for three using 
our broad-band-only SED fitting routines. Using Calzetti et 
al.'s (1997) nomenclature, these are clusters NGC 5253-1, 2 
and 6 (corresponding to, respectively, objects 129, 45 and 36 
in Table 2). For NGC 5253-1, Calzetti et al. (1997) derive 
an age of 8-12 Myr, while we find a best-fitting age of 14 ± 3 
Myr. Calzetti et al. (1997) and Tremonti et al. (2001) find 
ages of, respectively, 50-60 and 81q j Myr for NGC 5253-2, 
compared with our best age estimate of 5.0 ± 1.6 Myr. Our 
results thus support the younger age advocated by Tremonti 
et al. (2001), although our mass estimate, (6.4 ± 1.3) x 10 4 
Mq is somewhat higher than theirs, ~ 1 x 10 4 Mq. Note 
that although this cluster's derived mass tends towards the 
mass range where the effects of stochastic sampling become 
smaller, we caution that a small number of luminous red 
stars could easily affect its broad-band SED to mimic a 
younger age. Finally, for NGC 5253-6, Calzetti et al. (1997) 
find a best-fitting age of 10-17 Myr, which is a close match 
to our age estimate, 14 ± 2 Myr. 

We also compared our derived age and mass estimates 
with the equivalent results of Harris et al. (2004). Within 
the uncertainties associated with both sets of parameter de- 
terminations, our age and mass determinations of the dozen 
objects with parameter determinations in common are fully 
consistent with theirs (not shown). 

Recently, Harbeck et al. (2012) reported the discovery 
of three potentially massive clusters (M c \ > 10 5 Mq) with 
ages of order 1-2 Gyr in the galaxy's outer regions, based 
on HST/ACS imaging observations in the F415W, F555W 
and F814W filters. These clusters, which they refer to as 
#2, 6 and 7, form part of a more extensive population of 
intermediate-age and old clusters. From an initial focus on 
the galaxy's starburst core and its properties, more recent 
studies have attempted to explore the galaxy's overall star- 
formation history. The intermediate-age clusters of Harbeck 
et al. (2012), as well as our larger population of intermediate- 
age and old clusters support the notion that NGC 5253 is 
a very active dwarf starburst galaxy that has undergone 
multiple episodes of star and star cluster formation (see also 



Calzetti et al. 1997 and McQuinn et al. 2010a,b for similar 
conclusions based on studies of the field-star population). 

6.4 The cluster size distribution in the context of 
previous work 

Next, we consider the cluster size distribution shown in Fig. 
1 in the context of previous work. Harbeck et al. (2012) tab- 
ulate the FWHMs of the 28 objects in the galaxy's outer 
regions which they discuss (18 likely clusters and 10 back- 
ground galaxies). Their resulting size distribution (converted 
to Gaussian a's) is generally flat, however, and does not 
show a clear peak at smaller sizes similar to what we find 
in this paper. This is likely caused by their cluster selection 
procedure, which was based on visual inspection of the ACS 
images. In addition, we cannot directly compare our cluster 
parameters derived in this paper with Harbeck et al.'s (2012) 
results, because both samples are complementary, without 
any spatial overlap. 

Harris et al. (2004) also examined the NGC 5253 cluster 
sizes and include a figure showing their (FWHM) size dis- 
tribution as a function of cluster magnitude. Their size dis- 
tribution is based on analysis of their F547M image (which 
is closest in wavelength to the F555W filter we used). The 
vast majority of their sample clusters have sizes which are 
similar to the PSF size, with a small tail towards more ex- 
tended clusters. Although this trend is qualitiatively similar 
to our result, we point out that their size measurements were 
based on WFPC2/WF camera images, which have a factor 
of approximately 2 lower resolution than the ACS/HRC im- 
ages upon which we base our analysis. The bottom inset in 
Fig. 1 shows the equivalent figure for our data set, also us- 
ing the F547M magnitudes. The trend shown is similar to 
that found by Harris et al. (2004). We re-emphasize here 
that we specifically selected objects that appeared compact, 
which implies that we selected against including stellar as- 
sociations. 

Finally, we consider the cluster size distribution in a 
wider context. Bastian et al. (2005) examined the effective 
radius (_R e ff) distribution of their YSC sample in M51 and 
of the old GCs in the Milky Way. They represented the dis- 
tribution by a power law of the form N(R)dr oc r~ v dr, and 
found best-fitting power-law indices of r) = 2.2 ± 0.2 and 
2.4 ± 0.5 for the M51 and Milky Way clusters, respectively. 
This compares to rj = 3.4 for the YSCs in the merger rem- 
nant galaxy NGC 3256 (Ashman & Zepf 2001). Here, we find 
a best-fitting power-law distribution of the clusters' Gaus- 
sian <t's characterized by rj = 2.8±0.6, which agrees to within 
the uncertainties with these previous determinations - pro- 
vided that there is a one-to-one correlation between R e g and 
(Jg, which is a reasonable assumption if cluster profiles do 
not differ too significantly among galaxies. 



7 SUMMARY AND CONCLUSIONS 

Because of the wealth of available multi-wavelength data 
sets in the HST Data Archive, we chose NGC 5253 as target 
to compare two model approaches to determine fundamental 
star cluster parameters based on up-to-date input physics. 
We have gone significantly beyond previous studies of the 
NGC 5253 star cluster population. 
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First, although we still used medium-, broad-band and 
Ha imaging observations, we base our results on the highest- 
achievable spatial resolution. The gain in resolution com- 
pared to previous studies is of order a factor of two in both 
dimensions, while our accessible wavelength range tran- 
scends previous studies by incorporation of both the NUV 
and NIR passbands. Crucially, this allows us to at least par- 
tially break the age-met allicity and age-extinction degen- 
eracies. Second, we applied two spectral synthesis methods 
to our HST observations in up to 10 niters, taking into ac- 
count the updated effects of nebular emission and the im- 
proved physical understanding of the excess flux at red opti- 
cal and IR wavelengths found for a subset of young clusters. 

We used the most recent models to explore the origin 
of excess emission of a fraction of the YSCs in optical I and 
longer-wavelength NIR filters compared to 'standard' SSP 
models, assuming fully sampled stellar mass functions. Of 
the 182 objects in our sample with reliable photometry in 
at least 7 filters (i.e., uncertainties < 0.3 mag), we obtained 
best fits of sufficient quality (i.e., with sufficiently small xt 
and Q > 0.001) for 149 cluster candidates; we discarded the 
remaining 33 objects from our final sample. Of these 149 
clusters, 30 were characterized by a clear IR excess, while 
the remaining 119 objects exhibited SEDs as expected for 
fully sampled clusters. We compared our derived cluster ages 
and masses based on application of our novel Yggdrasil- 
based approach (Adamo et al. 2010a,b, 2011a,b) with the 
results from application of GALEV+ ANALYSED SED match- 
ing. The main features in the overall age and mass dis- 
tributions were well reproduced by both independent ap- 
proaches. Assuming a simple underlying Gaussian distribu- 
tion, the peak in the age distribution resulting from our 
GALev+AnalySED analysis occurs at (log(t yr -1 )) = 7.11, 
with o"(i ogt ) = 0.80. For comparison, the Yggdrasil-based 
analysis returned (log(t yr -1 )) = 7.10, with C(i og t> = 0.94. 
For the mass distributions, we find (log(M c i/M0)) = 3.38 
versus 3.37, and tX{i ogMcl ) = 0.46 versus 0.77, respectively. 

The NGC 5253 cluster population is dominated by a 
significant number of relatively low-mass (M c i <C 10 M©) 
objects. Assuming fully sampled stellar mass functions, they 
have ages ranging from a few x 10 6 to a few x 10 7 yr. This 
age range is in excellent agreement with the starburst age of 
the host galaxy (< 10 Myr). Most of our IR-excess clusters 
are found in the younger 'chimney' observed in the diagnos- 
tic age-mass plane based on our fully sampled SSP analysis, 
characterized by an age of ~ 5 x 10 6 yr and masses of up to a 
few x 10 4 Mq . The observational properties of these clusters 
are most likely affected by a combination of stochastic sam- 
pling effects and the presence of PMS stars. The IR-excess 
clusters in the older chimney, at log(t yr -1 ) = 7.1-7.2 are 
also of relatively low mass and, hence, we expect them to 
be affected by stochastic sampling effects, while at this age 
the stochastic colour variations owing to luminous red stars 
(including RSGs) will also become apparent. 

Our analysis also revealed the presence of a small num- 
ber of intermediate-age (~ 1 Gyr-old), ~ 10 5 M Q clusters, 
as well as up to a dozen old clusters resembling true GCs, 
in relation to both their ages (~ 10 Gyr) and masses ( > 10 5 
M©). Our detection of this fairly small number of true GCs 
is simply owing to selection effects and may not represent 
a genuine dearth of old, low-mass clusters in the galaxy. 
The presence of populations of young, intermediate-age and 



old clusters supports the notion that NGC 5253 is a very 
active dwarf starburst galaxy that has undergone multiple 
episodes of star cluster formation. The steeply increasing 
detection limit as a function of increasing age implies that 
even if there were large numbers of intermediate-age and 
old clusters of intermediate masses, we would not be able to 
detect them at the present time. 
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